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PREFACE 


An analytical technique for predicting the design requirements of the 
oombined eooustio end meohenloel loads lnourred during the Spaod Shuttle 
launoh phase dtie to the simultaneous oodurrenoe of low frequency transient 
inputs and aooustioally induced random vibrations has been developed* this 
dooument describes the underlying theOretloal rationale leading to the 
proposed design method* 
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1.0 INTRODUCTION 


During the laUnoh of a Spaoe Shuttle Tranaportation System (STS) the 
almultaheoua ooourrenoe of mechanically tranamitted low frequertey transient 
excitation, generated by the ignition Of the solid rocket boosters, and acous- 
tically induced random vibration requires that the design loads speoified 
STS payload components include the effect of this combined excitation. Docu- 
mented herein is a study that first develops a statistical description of how 
the transient-mechanical and acoustically induced loads are expected to 
combine, and then proceeds to describe and Validate a procedure for setting 
design requirements for loads expected to combine in this manner. A basic 
assumption in this study is that the individual acoustic load and mechanical 

load have been previously obtained. 

The combined load statistics are developed by taking the acoustically induced 
load to be a random population, assumed to be stationary. Bach element of 
this ensemble of acoustically induced loads is assumed to have the same power 
spectral density (PSD) , obtained previously from a random response analysis 
employing the given acoustic field in the STS cargo bay as a stationary random 
excitation. The mechanically induced load is treated as either (1) a known 
deterministic transient, or (2) a non-stationary random variable of known 
first and second statistical moments which vary with time. A method is then 
shown for determining the probability that the combined load would, at any 
time, have a value equal to or less than a oertain level. 

Having obtained a statistical representation of how the acoustic and mechani- 
cal loads are expected to oombine, an analytical approximation for defining 
design levels for these loads is presented using the First-Passage failure 
criterion. Using this criterion the probability of the first Occurrence of a 
oertain level (the design value for the odmbined load) within a prescribed 
period (the desired service life of the STS payload oomportent) oart be deter- 
mined. Empirical verification of this approach for establishing design loads 

is then provided. 
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While the procedure for establishing design loads is developed assuming a 
stationary aooustio pressure field, this study also Indicates how this pro- 
cedure can be modified to account for non-statlonary aooustio excitation. 
Finally, a method is presented for establishing design loads that acoount for 
random flight-to-f light variation in the aoouatically induced load. The 
salient features of a random process (the aoouatically induoed load) with 
whloh the reader should be familiar and that are employed in the text of this 
document are presented in Appendix A. 

2.0 A PROCEDURE FOR COMBINING ACOUSTICALLY AND MECHANICALLY INDUCED LOADS 

A procedure for predicting combined acoustically and mechanically Induced 
loads has been developed. The procedure assumes that the individual acoustic 
load and mechanical load have been previously obtained. TWo cases are 
considered: (1) the mechanically induced load is taken to be a known 

deterministic transient; (2) the mechanically Induced load is taken to be a 
non-stationary random variable. 

2.1 Deterministic Transient Mechanical Load Plus Stationary Random Acoustic 
Load 

In. this ease the mechanical load, X(t) due to a staging event, is taken to be 
a known deterministic transient (e.g. Figure 1). This information is obtained 
from techniques (computer codes) currently employed to determine the mechan- 
ically Induced transient load for staging events of expendable boosters, and 
entails a response Analysis of the pbyldad, as represented by a finite element 
model, to the known deterministic staging excitation. For any time, t^, as 
the.. response due to staging is taken to be deterministic, we have 


E 

K>] 

• x(t x ) 



E | 


« x 2 (t x ) 


(1) 

a 2 x(t x ) 

• * [ x S>] 

- E 2 [x(t x )J « 0 
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where 


e [ ] 


E 

E 


[«v] 

o 2 *^) 


indicates expected value, and nenoe 


la the mean value of X(t^) 


la the mean aquare value of X(t^) 
la the varlanoe of X(t^) 


Aa ahown In Appendix A, the aoouatioally Induoed load, Y(t), xa a random 
population aaaumed to be atationary and Qauaalan with: 


E 

E 


H 

[Atj] 


(Equation AW) 


o\ * E 


. / 

•'O 

[At)] - 


E 


S y (W) dW (Equation A18) (2) 

o 

[*/)] 8 E [y 2 U)] (Equation A19) 


where 


E 

E 


[ T(t >] 

[At>] 


S y (W) 


19 the mean value of Y(t) 
la the mean aquare of Y(t) 
la the Variance of Y(t) 

is the power speotral density of Y(t) (e.g., Figure 2) 


Explicit reference to time for the aooustio load statiatioal parameters will 
be eliminated in the sequel, i.e., 


E 


E 


[««] 

[f 2 (t)] 


w 

[y 2 ] , etc. 


Assuming a linear system, the oomblned load, Z(t), Is a Superposition (sum) 
of the loads due to the staging transient, X(t), and the aooustlos, Y(t). 


Z(t) 


X(t) ♦ Y(t) 


(3) 






i 

! 
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or at any given t. , t^ 
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X(t 4 ) 4 Y<t i ) 

taking expected values 

• 

E [Z(t i )] 8 

E [Z(t t )] 

X(tl) 0 

E [X(/ ± )] * E (/ 

x(t t ) 


That la, the mean value of the total load at any time, t ^ , la merely the 
load at that time due to the transient since the expected value of the load 
due to aooustics is zero. 


Squaring the total load 

z 2 (t i ) » x 2 (t i ) 4 Y 2 (t i ) 4 

and taking expected values 

x 2 < ti) 

E [z 2 (t i >] * e[x 2 ^)] + E 

where 

e fxdj) y ( t t )] * covcxap 

and COV (X(t i ) Ytt^) is the covariance of X(t A ) Y(t^) defined as: 

o 

COV(X(t 1 ) Y(t t )) s E {[XO^) -/E(X(t i ))] [Y(t t ) - B(Y)]} 
hence CCVW^) Y(t A )) • 0 

and, from Equation 7 E £z 2 (t A )J a X 2 (t i ) 4 cr 2 ( 

Now, the mean square value is the variance plus the square of the mean 
olitj 4 E 2 [Z(t t )] a X 2 (t 1 ) 4 q\ 


2 X(t. ' A ) 


X 

[j fi ] + 2E (xct|) Y(tj) 


X(ti) o 
Ydj)) + E [X^)j E \]t\ 


( 
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o|(t t ) ♦ X 8 (/ t ) • * °y 

* °X 


( 12 ) 

(13) 


•Ad the standard deviation of the combined load is equal to the Standard 
deviation of the aooustio load. 


tf 2(t 1 ) « (l 

Finally, as the acoustically induced load is normally distributed, we obtain 
the P** 1 percentile of the combined load from 


or, 


where 




4(t t ) * x u i> ’ ^ a r 


(15) 

( 16 ) 


X(t i ) is the meohanically induced load at time, t 1# due to the 
staging event, a ^ is the standard deviation of the stationary 
acoustically induced load and is obtained as the square root of 
the area under the aooustio load PSD determined in the random 
response analysis of the aooustio exoitation. K is the 
appropriate oortstant relating a multiple of the standard deviation 
to the P** 1 percentile for a normal distribution (e.g. for the 
95*** peroentile, K o 1.645). 

Thus, the predioted peroentile oombined load at any time, t^, is obtained 
by adding to the detemihistio load at time, t^, due to the transient* 

X(t^) , the appropriate oonstant multiplied by the standard deviation of the 
atatlonary aooustio load, 
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As presented above, the P th percentile level represents the probability 
that an ensemble member will, at a given instant, have a value equal to or 
less than this level, or equivalently were the phenomenon ergodio (which the 
oombined load is not), the proportion of time this member spends at or below 
this level. It should be noted, however, that this information about an 
ensemble member is not that desired for determining oombined load design 
values. Combined load design values are developed in seotion 3.0. The 
reason for including the present probablistic description for oombined loads 
is to support the design load rationale given in section 3.0. 


2.2 Non-Stationary Random Mechanical Load Plus Stationary Random Acoustic 
Load 


In the previous seotion we considered the problem of superposing the load due 
to acoustic excitation, whioh was taken to be a stationary random variable, 
on the load due to S staging transient which was taken to be deterministic. 

In this seotion we make the extension to the situation where the load due to 
the staging transient is not deterministic but is itself a non-stat i onary 
random variable (e.g. Figure 3). This situation arises, for example, when 
the predicted transient load comes from an ensemble of transient excitations. 

The losd, X(t), due to staging is taken to be a known, non-stationary random . 

variable. For any time, t t , we know £ X(t i ) and ^(t.) by methods 
currently employed in predicting the mechanically induced transient load for 
staging events of expendable boosters. This entails performing a series of 
deterministic response analyses employing, one at a time, the members of the 
staging transient excitation ensemble, and subsequently forming the 
Statistical param'/,ers of the resultant responses* — 

Note, as the transient load IS non-stationary 

g •* 


|*<vl 


a x(t 2 ) 


(It) 


J 

i--. 
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The results obtained in Appendix A pertaining to the load due to aooustios, 
Y(t), still hold - namely, for any tiftd. 


E [y] ■ 0 

E [y 2 ] 8 / Sy(W)dW U8) 

0 l 

°v ■ [•? *«-]• 

As before, the total load, 2(t), is a superposition (linearity assumption) of 
the loads due to the transient and the aooustios 


Z(t) * X(t) + Y(t ) 


(19) 


or, at any given time, 
ZCt^ a x(t t ) ♦ 
taking expected values 
E [ttt,)] * E 

E [Z(t 1 )3 a E 


Y(t t ) 

o 

tX(t A )l + E \/f\ 
[X(t i )] 


( 20 ) 


( 21 ) 


That is, the mean value of the total load at any time, t^, is equal to the 
mean value of the transient load at time, t^, since the expeoted value of 
the load due to acoustics is zero. 


Now, squaring the total load and taking expeoted values - 


E [Z 2 (t 1 )] « E tx 2 (t i )] ♦ 

where 

E [x(t 1 )Y(t i )] a COV (X(t A ) 

and the oovariance is: 

COV (X(t t ) Y(t i )) a E { t X(t ± ) 


E [Y 2 ] a 2E [X(t 1 ) Y(t A )] (22) 

o 

Ytt^) + E [X(t 1 )] E \jf\ (23) 

- E (X(t A ))] [Y(t i ) - E(Y)]} (2k) 
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which can be written in terms of the correlation coefficient, Px y as* 

cov mt^ m^) ■ 9 x(t t ) a t 

where - 1 * P *<t> *<*> * * 1 

hence., equation 22 becomes 

■ [**(»,)] . ■ [«**,>] * E n * * p «v«v ° x< v * 

— tauare value is the varianoe plus the square of the jwan 


again, the mean 

4 (t t > * 


h i,uafe ealu. le th. eerlano. piu* -r~ - o 

* * 2 [*(<)] • ♦ e2 M ‘ °« 4 E $ <27 ’ 

oancela 


♦ 2 P 


X(t 1 )Y(t i ) ff X(t t > Y 


®« tl , 8 ‘5tt 1 > + ** * 2 P *<V «V tfx( V * Y 

and the standard deviation of the combined load is: 

• Wl > ■ [<»«*,> ♦ »S « 2 p «* t > «v ° x( v ° X J <29) 

nnall „ a, before, a.,u.ln« the tote! load U non»U» bletrlbeteb. the 
p th peroentile is given as: 

- p t P* <30 

■ E [ zct i'] * * 

ot V 

p* Ui> . t * * r * ♦ 4 * 2 p xtt 1 mt 1 ) "xftj) °jJ (W 

To hake further progr... .. «»t «... ao». ..auction re,ardl«« th.je.re. of 
correlation between the random variables X(t) and Ytt). e oons 


oases i 


1. If the load due to th. treneleht 1. totell, unodrreleted 
with th. load due to aoouetloi, i.e., X(t)Kt) • 0 
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then 

% ! 

*’<*!> • e [ x(t l>] » ** [*«,,) * a l] <3*> 

That is, the standard deviation of the total load la a 
Root-Sum-Square (RSS) of the standard deviations of the 
transient load and the acoustic load, 

2, if the load due to the transient is perfectly correlated, 
with positive correlation, With the load due to acoustics, 
i*e., ^X(t)Y(t) s+1 (the loads are in-phase and 
reinforce each other) 

then 

P Z(t 1 ) • E [*<v] ♦ t'* < °x( tl ) ♦ <V (33) 

That is, the standard deviation of the total load is the 
sum of the standard deviations due to the transient and 
the acoustic loads. 

3. If the load due to the transient is perfectly correlated, 
with negative correlation, with the load due to acoustics, 
i.e., ^x(t)Y(t) 8 *1 (Tne loads are out of phase and 
subtrdot from each, other) 

then, 

P 2 (t 4 ) * E [ X( V] 4 K * | a x(t 1 ) * *y| (3<0 

That is, the standard deviation of the total load is the 
absolute valiie of the difference of the standard deviations 
due to the transient and the acoustics. 
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As the aooustioaliy and meohanleally induced loads at Shuttle lau*'^' arise 
from fundamentally two different exoitatlon souroes and arrive at the payload 
by different and olroultous paths, it seems reasonable to take the first ease 
( PXY ■ 0) as being representative of the load for that event* 

3*0 COMBINED DESIGN LOADS FOR FIRST PASSAGE FAILURE 

In the previous sections it has been demonstrated that, under the assumptions 
made, the combined load oan be represented as a non-statlonary Cauasian random 
variable with time-varying mean* As given in equation 5 the time-varying mean 
of the combined load Is merely the time-varying deterministic load due~to the 
mechanical transient. As given in equation 14 the standard deviation of the 
combined load Is the standard deviation of the stationary aooustioaliy induced 
load. The probability that a member of the combined load ensemble would, at a 
given instant, have a value equal to or less than a oertaln level, was also 
determined. It was noted that for an ergodio phenomenon this oould be 
considered as the proportion of time this member spent at or less than a 
certain level. It might be expected that a combined load design level oould 
be considered adequate if set such that this proportion of time was large. 
However, the proportion of time spent at or less than some level is not 
particularly relevant slnoe a failure could ooour the first time the 
magnitude of the combined loads equaled, or slightly exceeded, the design 
level. Therefore, for design purposes, one is Interested in the first time 
that a oertaln level is reaohed and what the probability is that this will 
ooour during the service life of, for example, a payload. This time will be 
different for each member of the oomblned load ensemble. The ensemble of 
M first M times has some distribution whioh must be known to establish the 
desired probability of ooourrenee of a given value for the oomblned load. 

The above considerations fall Under the oategory of "First-Passage" theory. 

A classical problem in stoohastio theory is to determine the probability P(T) 
dtP that the value of a random prooess surpasses a. threshold for the first time 
during the Interval from T to T ♦ dT. The resulting First-Passage probability 
density P(f) has considerable importanoe as a reliability measure in random 
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vibration studies. However, an exaot solution to this problem has not been 
found for even the simplest version of the problem - whioh is to oonslder the 
stationary response of a simple oaolllator exolted by white noise. However, 
a number of papers dealing with approximation methods whioh provide 
First-Passage probabilities have been written. 

In the study of many random prooesses the mean is assumed to be aero. If the 
random prooess serves as input to a linear system, a non-zero or time-'rarylhg 
mean oan be oonsldered separately and handled as a- -deterministic prooess. 
Thus, as seen in the previous sections, if one is only Interested in response 
statistios, such as mean and standard deviation, a time-varying mean presents 
no problem. However, the First-Passage problem for a system subjeoted to 
random exeitatlon is not merely one of finding the response statistios. 

In general, the first-passage probability density depends in a oomplloated 
manner, on the characteristics of the dynamlo system involved, on the nature 
of the exoitation, on the initial conditions Imposed, as well as on the 
magnitude of the ..threshold* 

The first-passage problem of interest at hand is that for a structure sub- 
jeoted to random exoitation with time-varying mean. That is, it has been 
shown that the combined load oan be taken as a deterministic prooess, the 
mean (l.e. the mechanically induced load), plus an additive random prooess 
with zero mean (the aooustlcally induced load). Having previously obtained 
the statistios for the oombined load, a method for approximating the first- 
passage probability is now presented. 

The time-varying mean is handled as a deterministic prooess and allows for 
formulation of the first-passage problem for a fixed barrier from two differ- 
ent outlooks. The first outlook is the one of considering tho oombined load 
as a random prooess with time-varying mean and studying first-passage statis- 
tics for a fixed barrier (e.g., Figure 4). The seoond view is to oonslder 
the oombined load as a random prooess with zero mean and to Investigate the 
first-passage problem for a time-varying barrier (e.g. Figure 5). In this 
oase the barrier is. made time-varying by subtracting the prooess mean. The 
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* 

two outlooks are In essence equivalent but the latter Is adopted here* Thus, 
the first-passage problem is considered in terms of a time-varying barrier 
for a aero mean stationary Gaussian random process. 

The simplest approximation to the failure rate, ealled the Poisson approxima- 
tion , assumes the barrier ero.. '.ngs ooour so rarely that they oan be consid- 
ered as statistically Independent events. The assumption of independence of 
crossings has been attacked but it has been pointed out that any tendertoy of 
the crossings to cluster in groups or clumps (l.e*, not be statistically 
independent) makes the Poisson assumption conservative, and hence the use of 
the Poisson assumption in design introduces an error on the conservative side. 


The Poisson process is one whose properties have been thoroughly explored. 

In particular, it is known that for Poisson processes the First-Passage prob- 
ability density function is an exponential function with a single parameter 
which is simply related to the expected rate of crossings, v & . That is: 



Thus, on the basis of the assumption that the crossings of the level Z(t) « A 
constitute a Poisson process we have a "solution” to the distribution to the 
time of "allure. When the distribution density, P(T), for the time to failure 
is known, the probability of failure in the interval T^ < T < Tg is: 


Pttg,^) 

or, 

PtTg.Tj) 



PCT) dT 


exp 



(36) 

(37) 


For a positive barrier, a(t)>0, the failure rate is the upcrosSlng rate, 

V + , of the level a(t) . For a negative barrier a(t) < 0, the failure 
rate is the downcrossing rate, V & ”, of the level ot(t). In general the 
total expected rate of barrier crossings is the sum of the upcrosslng rate 
for a positive barrier and the downorossing rate for a negative barrier, i.e.: 
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The barrier crossing rate, i> a , for a atationary random prooeaa for a 
fixed barrier was firat derived by Riee (Referenoe 1). Thia waa extended to 
the oaae of general ourve oroaaing of a non-atationary Oauaaiari prooeaa by 
Cramer and Leadbetter (Referenoe 2). The uporoaaing rite, v*(t ,<*)» for a 
poaitive barrier level and the dovmoroaaing rate, v"(t,tt) for a negative 
barrier level are given by (Reference 2, page 288). 

* -^7='* ,xp + n * <r,) l 


(38) 


v?t,a) * 


. 0(1 - uV 


•J Ta 


exp 


y) bt? exp w ■ 


n(l - *(n)) 


where (for the problem at hand): 

8 B C^(t)] ; tHe variance of the acoustioally induced load 

£ 2 8 E Ci 2 (t)J ; the variance of the time rate of change of the 

acoustically induced load 

M « E r?(t) Y(t)1 : The normalized covariance or correlation 

a(t) a(t) coefficient for the above two parameters 

a(t) .is the time-varying barrier, at(t) * A - X(t) (i.e. o(t) is 
obtained by subtracting the meohanically induced load from the 
fixed barrier, A). 


Further, we define: 


i) « 


a ■ « 


(• 


(1- M ) 1/2 

4 - (° (t) ) 


dt 


a 


* j p) 


« - X ( t ) 
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<f>(rj) designates the normalized Gaussian distribution function: 


4>(TJ) 


yi # - 




For the problem at hand, the aooustioally lnduoed load is taken aa 
stationary. It is known (Reference 3, page 33) that for a stationary 
prooess, Y(t) there is no correlation between Y(t) and Y(t), i.e., 

E [Y(t) Y(t)] a 0 


and hence, 


E [ Y( t ) Y(t ) 

a(t) STt) 


Further, for a stationary process, 


.2 

a * 


’Sy(W)dW 


Thus, for the general curve crossing of a stationary Gaussian process, 
equations 38 become: ^ 

K + Ct,a) = jt, exp (^7) f estp (jps) ' I / e' tS,2 dt 

•00 

* & § «*P (^5) [ exp (^) ■ ? / e ' ,2/2 dt * 'fi* f 


Further, should the stationary Gaussian random prooess with zero meah also be 
narrow band - with center frequency W Q then it is known (Reference 3 » page 
44) that:.. 


s _a 
a 
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Hence • the curve crossing of a stationary, narrow band, Oausoinn process la: 

/(t ,«) 


exp 


of a stationary, nsr^ 

(i) ; xp (i§?) ' 


-m 0 o _ A 

wr » / 


& 

T l 

V »l 


0 *' 12 « 


.« ■ % exp ($) [ exp (s£i) '«>£ *’* * * + ^ CT '“o ] 


(40) 


Finally, should the barrier be constant, 


o(t). * a; a a 0 


Then the terrier crowing of . fixed terrier by e et.tlon.ry, nerrow b.nO, 
Gaussian process is given by: 

v*(a) * iT(a) * Us ex P (=^£\ (“1 

2lt \ I 


which is the result first given by Rloe (Reference 1) or (Reference 3, P»g« 
107). 

For the First-Passage problem of the combined load, under the assumptions 
stated t equations 39 are in general applicable. If in addition the 
acoustically induced load can be taken as narrow band, then equations 40 may 
be used. In either case, equations 39 or equations 40 are employed to 
calculate the expected rate of barrier crossing, V & 8 v+ ♦ V~ , which 
is substituted into equation 37 to evaluate the relationship between the 
probability of failure P(T 2 ,T 1 ) during the interval Tj<T<T 2 and the 
barrier level A (where A is the combined design load). 

To illustrate, consider the somewhat restrictive case in which the acousti- 
cally induced load can be taken as a narrow band process of center frequency 
W and standard deviation a. This, and other, information on the acousti- 
cally induced load is obtained from a random response analysis with the 
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structure represented by a finite element model and employing tre acoustic 
pressure field as a stationary atoohastio excitation. Further, the meohani- 
oally induced deterministic transient load, x(t), is previously obtained by a 
response analysis of the payload, as represented by a finite element model, 
to the known deterministic transient staging excitation. Assume the oombined 
loading oondition lasts for 10 seoonds (l.e., « 0 and T g » 10). We 

seek the barrier level, A, such that the probability of failure, P(T 0 ), 
during this interval is one percent (i.e., the probability of suooess is 99 
peroent, P auoeess * 100-P fai ^ ure ). 


The relevant equations are equations 3/ and 40: 



where T = 10, P(TJ a 0.01 

o o 

o(t) - A - X(t), a(t) a - i(t) 

and W Q , X(t) , X(t) , a are known a priori and we solve for A (the oombined 
design load). 

An interpretation of the above is that out of 100 flights, on the average, 
only one will exceed the design level A when exposed to the combined load of 
ten-second duration. 



Note that intuition ia aerved in that the above equations for predicting 
flrat-padOage failure are funotione of amplitude* duration, and-frequenoy 
oontent of the oombined loading oondition. The longer the-duration of the 
oombined loading oondition the greater the probability that a preacribed 
barrier level, A, will be exoeeded. The greater the frequenoy oontent of the 
oombined lo ad i ng oondition, the greater the number— of excursions per unit 
time and henee the greater the probability that a presoribed barrier level 
will be exoeeded. 

Obviously, the solution of these equations must be carried out numerloally 
employing a digital computer. The problem with any type of orossing rate 
approach for a rapidly moving barrier is that the statistics must be computed 
at frequent intervals. This can beoome quite expensive, especially when the 
frequenoy is high. 

Extension of this technique for First-Passage failure to the case where the 
meohanloal load is taken as a non-stationary random variable is made by 
fr. piHng the oombined load mean as given in equation 21 and the combined load 
atandard deviation as given in equation 29 and employing equations 3_8 for the 
orossing rates. 
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4.0 EXPERIMENTAL INVESTIGATION 


To evaluate the proposed method of combining low frequency transient loading 
with aooustloally induced- loads, experimental methods were employed whereby a 
test specimen was separately and simultaneously exposed to aooustlc 
excitation and portable shaker force transients. 

The test speoimen (Figure 6) was a 14-foot diameter right-circular oylinder, 

10 feet long, with double wall wooden end caps. The cylinder was of .125 
inch aluminum akin with circumferential and longitudinal stiffeners, and was 
tested with a simulated payload component that weighed 45 lbs and had an 
installation first mode at 27 Hz (Figure 7). The speoimen was suspended by 
steel cables from the work platform in the large reverberant chamber of the 
LMSC acoustic test facility. A shaker was suspended using steel cables and 
the stinger attached to the mass simulator, as shown in Figures 6 and 7. 

The test instrumentation was comprised of: (1) four cell microphones to 

monitor the acoustic field present in the reverberant chamber; (2) twelve 
Endevco 2220 accelerometers to monitor the mass simulator responses located 
as shown in Figure 8 and 9{ and (3) four uniaxial strain gages (two pair of 
back-to-back gages) as shown in Figures 8 and 9. In addition, a load cell 
(Figure 7) was employed to monitor the shaker input force to the specimen. 

It should be noted that the techniques discussed herein are methods for 
combining .signals and subsequently setting "design levels" for those signals 
under a first-passage failure criterion, regardless of the physloal signifi- 
cance of the signal. Thus, while the variable of immediate interest is design 
load, or stress, for the purpose of empirically verifying the proposed tech- 
niques more readily measureable variables such as acceleration or strain, will 
serve. 
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Test Configuration - Close-up 





rrar* 










AS the first-passage failure level is a function of the amplitude! duration 
and frequency content of the oombined response! two aooustlo spectra and 
three shaker transients were employed in the testing. The two aooustlo 
speotra! differing in low frequenoy energy oontent. are shown in Table I. 
Based on the results of a shaker sine-sweep test! the shaker transients were 
seleoted that would engender strong mass Simulator responses 6f magnitudes 
oomparable to the responses generated by the aooustlo speotra. The transient 
exoltatlon applied by the shaker was oontrolled "open-loop" using a magnetic 
tape loop eontainlng the desired transient. One tape loop was made 
containing the first transient whioh was repeated 100 times at 10-seoond 
intervals. The mathematioally constructed transient was of 4-seoond duration 
with a sine envelope and frequenoy oontent near a resonance (27 Hz) of the 
mass simulator to ensure suffiolent response amplitudes. Two additional tape 
loops were made containing the second and third transients! eaoh also 
repeated 100 times. The second transient, repeated at 10-se0ond Intervals 
was of 4-se6ond duration with a Sine envelope and frequenoy of 50 Hz. The 
third transient, repeated at 2.5-second intervals, was of 1-seoond duration 
with a sine envelope and frequenoy at 270 Hz. 

To confirm repeatability of the shaker transient excitation and to obtain the 
response. X(t), due to the meohanioal transient, eaoh transient in turn was 
applied 100 times and the response of the twelve aooelerometSrs and four 
strain gages recorded. To obtain the response ¥(t), due to the aooustlo 
fields, eaoh speotrum (Table I) was applied to the speoimen in the 
reverberant chamber for a duration of three minutes and the accelerometer and 
strain gage data recorded. Finally, to obtain oombined responses, 2(t). 
three combinations of aooustlo field and shaker transient, as listed below, 
were tested and the aooelerometer and strain gage data reoorded: 

1. speotrum A plus transient 1 applied 100 times for a duration of 1000 
seoonds 

2. speotrum B plus transient 2 applied 100 times for a duration of 1000 
seoonds 

3* speotrum B plus transient 3 applied 100 times for a duration of 250 
seoonds 


6431a (04/13/83) 
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TABLE I 

ACOUSTIC TEST SPECTRA 





ORIGINAL PAGE IS 
OF POOR QUALITY 

Preliminary data processing was completed for one accelerometer and on# strain 
gage for the teat sequence. The aeotleroneter aeleoted was acoeierdaeter #6, 
as it was a normal accelerometer located in the corner bf the mass simulator 
(see Figure 9), and hence had a strong response signal for both simulator 
translational modes (e.g. 27 Hz) and rotational modes (e.g.i 50 Hz). Due to 
geometric symmetry, the strain gage selection was arbitrary and gage #2 was 
chosen. 

For the transient-only tests employing transient #1 (27 Hz) and transient 02 
(50 Hz) both of which were of 4-seoond pulse duration, the response data X(t), 
were digitized at a rate of 2,000 samples per seoond, or 8,000 samples per 
pulse. For the transient only test employing transient #3 (270 Hz) which was 
of 1 -seoond puls© duration, the response data were digitized at 4,000 samples 
per second (pulse). The time derivative of these response data were then 
calculated using a simple "delta-delta" formulation, i.e., 

X * AX/ At 

For the short duration high frequency thira cansieirt, the strain gage data 
were low in magnitude and of questionable quality. Hence, for the high 
frequency (270 Hz) transient, only the accelerometer data was processed. 

For the acoustic only tests, the power spectral densities of the responses, 

S (W), were determined. From these, the standard deviations of the 

y 

acoustioally induced response, <jy* and the time rate of change of the 
acoustically induced response a», were calculated 



For the combined acoustic/transient tests employing transient #1 and transient 
02, response data, Z(t), were digitized at 2,000 samples per second (8,000 
samples per pulse). For the combined aooustic/transient test employing trans- 
ient #3t the accelerometer response data were digitized at 4,000 samples per 
seoond (pulse). 
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Pop the adoustio-only testa and the combined aooustio/ transient tests, coll 
microphone one-third ootave SPLa were generated to ensure repeatability of 
the aooustio fields* 

H.l Comparison of Theory and Experiment 

Comparison of the first-passage level theoretical prediction and experimental 
data for the oomblned aooustio/ transient loading conditions are presented in 
Figures 10-14 and Tables II - VI as described below* Employing a digital 
computer, the theoretical prediction computations were performed following 
the general steps 2 through 5 of Table VII. The empirical curves were 
obtained by merely observing the percentage (number) of a response ensemble 
that did not exceed (even one time) each prescribed barrier level* 


Figure 

Table 

Instrument 

Condition 

10 

II 

Aocel #6 

Spectrum A + Transient 1 

11 

III 

Accel #6 

Spectrum B + Transient 2 

12 

IV 

Aocel #6 

Spectrum B + Transient 3 

13 

V 

Gage #2 

Speotrum A + Transient 1 

14 

VI 

Gage 92 

Spectrum B + Transient 2 


The theoretical predictions are made employing equations 37 and 39: 

P(T) = J v a exp ( -y* r a dr) dt 


where, v a = 


v + + v " 
a a 


and 


a 


--h 

' -00 


' -h i ex P 




(37) 


(39) 


is the variance of the acoustically induced response 


a 2 


is the varianoe of the time rate of change of the 
acoustioally induced response 
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ii the timc-Varying barrier, n(t) ' * - x(t > 
o (t) la obtained by aubtraotlng the moehanioally 
Induced reaponae from the fixed barrier) 


« ■ d/dt ( ft (t) ) » -x 

t • H or 1 (i.s. the eombined load la. applied for either 
it second a or 1 aeoond) 


Prom th. transient-only run. are obtained th. transient response. *<»>, 

Z tlm. rate of change of th. tr.n.l.nt response, X<t>. Prom th. ao«..tlc- 
Z, Z* Z Obtain* C . th. .tanl.ro Elation of th, acoustically lhhuo.0 

reaponae and 6. th. .tanddrd deviation of th. tlm. rat. of the 

acoustleelly lnduoed respon... Th... tuo ..... permit th. Predicted rin 
passage probability of failure to be made for various levels, », * 0 •*“ 

37 and 3W. The probability of not crossing a Riven level 1. the eompllmont 

of th. probability of failure given by eouatlon 37, i.e.. 


'SUCCESS * 100 • 

To the ..tent that th. proposed approaoh Is valid, the eombined response, 
ulll follow the first-passage probability given b, equation, 37 «>d 3.. 

„ seen In figures 10-1« and Tables Il-Vl, the agreement between 
paasags failure theoretloal prediotlon and th. ^lfio.1 data is quit, g . 
With th. ..ception of th. strain gag. data for the speetrum ' 

eort.ln.tlon, th. predictions are .lightly oonservat.'vs showing typical posl- 

- n M i t this result is expected and is attributed to 

tive margins of 1.0 to 1.5 dP. This result > . 

th. pr.vlou.ly stated conservative Poisson assumption that the barrier 

l«s are statistically independent events. Th. strain gage data for th 

apeotrum B/tr.n.l.nt * combination (Plgur. U and T-bl. VI) 

agreement between theoretical prediction and test data but without 

the difference between prediction and test balng at most b-'bB.Whynom, . 

ein i. paroelved for this o... ha. not be.n axplalnad at this 

th .„ to point out that the dlaparltla. between prediction and « plrl0,l d * 

are -.11 within th. historical unoertalntla. of typical vlbrcacoustl. data. 
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Fig. 10 First Passage Level Accelerometer 6 - Spectrum 
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Pig. 12 First Passage Level Accelerometer - Spectrum B/Transient 3 













TABLE X! 

FIRST-PASSAGE LEVEL COMPARISON 
Accelerometer #6 
Spectrum A ♦ Transient 1 
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TABLE III 

FIRST-PASSAOE LEVEL-COMPARISON- 
Accelerometer #6 
Speotrum B 4 Transient 2 


Probability of 
not Exceeding 
(*) 

Theoretical 

Prediction 

<g‘s) 

Empirical 

Data 

(*’8) 

Margin 

Ratio 

<dB) 

10 

19.2 

16.5 

1.16 (1.3) 

" 20 

20 .0 

17.3 

1.16 (1.3) 

30 

20.5 

17.6 

1.16 (1.3) 

MO 

21.0 

17.9 

1.17 (l.M) 

50 

21.5 

18.5 

1.16 (1.3) 

60 

22.1 

19.1 

1.16 (1.3) 

70 

22.7 

19.9 

1.1M (1.1) 

80 

23.5 

20.8 

1.13 (1.1) 

90 

2M.7 

22.2 

1.11 (0.9) 

92 

25.0 

22.2 

1.13 (1.0) 

9M 

25.5 

23.1 

1.10 (0.9) 

95 

25.7 

23.1 

1.11 (0.9) 

97 

26.5 

2M.2 

1.10 (0.8) 

99 

27.9 

2M.7 

1.13 (l.D 
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Probability of 
not Exoeeding 

ill 

Theoretical 

Prediction 

10**® inohe a/inch 

Emplrioal 

Data 

•6 

10 4 nehea/inoh 

10 

1071 

865 

20 

1111 

916 

30 

1190 

996 

90 

1166 

979 

50 

1191 

1011 

60 

1219 

1053 

70 

1298 

1090 

80 

1287 

1128 

90 

1393 

1192 

92 

1361 

1193 

99 

1383 

1203 

95 

1399 

1208 

97 

1931 

1221 

99 

1998 

1290 


I 


Margin 

Ratio 

(dB) 

1.29 (1.9) 
1.21 ( 1 , 7 ) 
1.21 ( 1 . 6 ) 

1.19 (1.5) 

1.18 (1.9) 
1.16 (1.3) 

1.19 (1.2) 

1.19 (1.1) 
1.13 (1.0) 

1.19 (l.l) 
1.15 (1.2) 
1.15 (1.2) 
1.17 (1.9) 
1.21 (1.6) 
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TABLE VI 

FIRST PASSAGE LEVEL COMPARISON 
Strain Qage #2 
Spectrum B 4 Transient 2 
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Probability of 
Not EXoeeding 

Theoretloal 

Prediction 

Empirioal 

Data 

KIM 

(%) 

-6 

10 lnohes/inoh 

10“ 6 lnohes/inoh 

(dB) 

10 

A 61 

A5A 

1.02 

(0.13) 

20 

A78 

A69 

1.02 

(0.17) 

30 

A95 

A90 

1.01 

(0.09) 

AO 

509 

500 

1.02 

(0.15) 

50 

521 

513 

1.02 

(0.13) 

60 

533 

528 

1.01 

(0.08) 

10 

5A6 

SAO 

1.01 

(0.10) 

80 

567 

569 

1.00 

(-0.03) 

90 

592 

598 

0.99 

O 

• 

O 

1 

92 

597 

611 

0.98 

(-0.20) 

9M 

609 

623 

0.98 

o 

CM 

• 

O 

9 

95 

618 

62A 

0.99 

(-0.08) 

97 

636 

637 

1.00 

(-0.01) 

99 

6A5 

6A8 

1.00 

(-0.OA) 
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Having achieved rather good agreement between the theoretloal first-passage 
prediction and the empirioal data, it ia not aurprialng that the aaaumptlona 
upon whioh that prediotion ia baaed-4iere reasonably valid for the teat 
configuration. A brief dlaouaaion of the empirioal verification of the 
underlying aaaumptlona la found in Appendix B. 

5.0 AN APPROXIMATION TO A NON-STATIOMABY ACOUSTIC PRESSURE PIELD 

The prevloualy developed analytical technique for predicting the combined 
loada during STB launoh due to tranalent meohanloal inputs and acoustically 
induced random vibrations ia baaed upon three assumptions: (1) that the sys- 
tem (struoture) is linear thereby permitting superposition of the two loads; ] 

(2) that the aooustio load IS Gaussian; and (3) that the aooustio load is j 

stationary. During STS launoh the third assumption of the aooustio load 
being stationary Is not physioally justifiable, but as the true evolutionary 
nature of the aooustio pressure may not be known, and in faot the aooustio 
forcing function often is. taken as a stationary envelope of the evolutionary 
aooustio SPL, the assumption of the aooustio load being stationary is made as 
a. pragmatic analytical expedient. 

Should the evolutionary nature of the aooustio pressure be known, this truly 
tton-statlonary forcing funotion could be approximated as being quasl- 
stationary in temporally contiguous intervals. The previously developed 

analytical techniques of oomblning transient meohanloal inputs with aooustl- 

oally induced random vibrations and subsequently predloting the probability 
of flrat-paaaage barrier oroasing of a prescribed level, A, oould then be 
affeoted for eaoh individual lntorval. The probability of flrat-paaaage bar- 
rier oroasing of a level, A, for the entire STS launoh event is then obtained 
as the sum of the probabilities of flrat-paaaage barrier oroasing of that 
level, A, in eaoh of the sub-intervals comprising the launoh event. This 
quasl-stationary approximation la illustrated in Figure 15. 







42 


RMS PRESSURE 





ORIGINAL PAGE 18 
OF POOR QUALITY 


6.0 THE! EFFECT OF FLIQHT-TO- FLIGHT VARIATIONS IN THE ACOUSTIC PBESSUBE FIELD 
PSD ON THE PREDICTED COMBINED LOADING 


The question may anise as to how flight-to-fiight variations in the acoustio 
field power speotral density , PSD, should be handled. To establish an 
aoouatic load based on the average of the flight-to-f light variations appears 
to be a non-conserVative approach but to base it on the maximum seems overly 
conservative. This seotion proposes an approach to handling this variation. 
Consider an ensemble of PSD* a representing the flight acoustic pressure 
fields (e.g. Figure 16). For a spatially uncorrelated aooutlc pressure field 
(i.e., negleoting the crosspower spectral density Of the acoustic pressure), 
the response power spectral density is related to the aooustio pressure power 
spectral density by the transfer function, T^ , J 

S y (X lf W) * T tJ Sp<Xj,W) (*»2 

where, T^ « Htt^X^W) H«(X it Xj,W) 
and, 


s y (x lf w) 


is the PSD of the response variable y, at 
frequency W and location X^. 


H(x it XV) is the frequency response function relating 

the response at location X^ to an acoustic 
pressure at location X^ at frequency , W. 


H«( x i, Xj,W) is the complex conjugate of H(X it X^, W) 

Sp(Xj,W) is the PSD of the acoustic pressure at 

location Xj and frequency W. 


Assuming the field is stationary, (i.e. statistical parameters are independent 
of time) and taking expected values across the ensemble of flight acoustic 
pressure fields results in, from equation b2: 


1 

i 


4M 





l 












m 


E 
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Squaring equation 42 and taking expeoted values across the ensemble 
results in: 


( 43 ) 



(44) 


As the mean-square of a random variable is equal to the varianoe plus 
the square of the mean, equation 44 may be rewritten as: 


r S v (X i ,W) 


+ E 


SJXj.W) 
* 






o 2 S p (X r W) 


+ T, 


ij 


E 2 


S p (Xj,W) 


(45) 


which by equation 43 reduoes to 


2 _ 2 2 
° S y (X if W) 8 T ij a S p (Xj,W) 

and taking the square-root of equation 46 

< S y (X i ,W) 8 T iJ °S p (X jt WJ 


( 46 ) 


( 47 ) 


where 


(X ,W) and °S (X ,W) are the fli Sht-to-flight standard 

deviations of the PSDs of the response variable, Y, and 
the acoustic pressure, respectively. 


Equations 43 and 47 state that the expected value (flight-to-flight mean) of 
the response PSD and the flight-to-flight standard deviation of the response 
PSD are obtained from mapping the expected value of the acoustic pressure PSD 
and the standard deviation of the aooustlo pressure PSD, respectively, 
employing the transfer funotlon, 



» H(X A ,XyJW) H«(X 1 ,X J ,W) 
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Now assume that the c^pui^jlgr^of acoustic pressure PSDs (from flignt to 
flight) are normally distributed. If the aooustie pressure PSD is normally 
distributed, then so is the response PSD, for the latter is linearly related 
to the former (equation 42). With the ensemble of response PSDs normally 
distributed, the poroentile of the response PSD is given byt 


<(,,, • E [VV»] * 

where k M * is the appropriate constant relating a multiple of the 

f u 

standard deviation to the M percentile for a normal 
distribution. 


(48) 


Substituting equations 43 and 47 into equation 48 results in? 


<«,« * T U E [VV°] * ^ [ T iJ 


or 


(49) 


P M 

p s y <* 1 .'0 


8 T 


ij 


E [s p (Xj,W)] 


♦ k” 


S p (X j ,W)| 


But i with the acoustic pressure PSDs normally distributed! the term in 
brackets, is the M th percentile of the acoustic pressure PSD. Thus, 



.V' 


M* 

8 T ij P S p (Xj,W) 


(50) 


Equation 50 states that the M th percentile PSD of the response variable is 
obtained by mapping the M th percentile PSD of the aooustie pressure field 
by the usual transfer function, * IKX^X^W) H*(X^,Xj,W). 
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The standard deviation of the stationary random acoustical)./ induced load, 

* 0y, ia obtained aa the square root of the area under the acouatio load PSD 
determined in the random reaonae analyaia of the acoustic excitation, 


ff V(X,) 



8 y (X,,W) 



(51) 


M 

When the standard deviation of the acouatioally induced load, a”(X 4 ), 
ft X X 

la based on the M percentile of the acoustic pressure PSD, equation 51 
becomes: 


a 


M 

Y(Xj) 



M % 

S y (X it W) 


1 


dW 




or from equation 50 


M _ T /•* p M % 

9 you J A ij p s D (X r w) 

o p 1 


-ll 



(52) 


The choice of which percentile, M, to use is a subjective, non-analytical 
decision. 


Recall the case when the mecnanically induced load, X(t),. was tateen-to be a 
known deterministic transient.. For this case the P th percentile of the 
total load, Z(t), was given by; 


P Z(t) 8 X(t) + (53) 

Now, when this percentile of the total load is in turn based upon the 

th m 

M percentile PSD of the aooustic pressure, substituting <jy for <Ty 

(with the location 1 dependence implied) equation 53 becomes: 


P 


% 

Z(t) 


X(t) + K 1 


M 


(54) 


48 


ORIGINAL PAGg 13 
OH POOR quality 


or 



X(t) + 



,M % 

yw> 



(5A ) 


Similarly, for the case where the mechanically induced load la a non- 
stationary random variable, the P th percentile of the total response, 
Z(t A ), at any time, t if was given by: 


p l<v - B 


|*X(tj)l + K P o^L. v + <r| + 2p _ 1 

L *J [*V Y P X(tj) YCtp ‘Xttj) a Yj 


And, when this p percentile of the total response is in turn based upon 
the M percentile PSD of the acoustic pressure, (i.e., replacing a 
with oh Y 

1 


zctp 


E 


[*v] + kP% 4«, 




*li 

} + <r*Y + 2p X(t i >Y(t i ) °X<t i ) °Y 


(55) 


M 

where in equation 55, Cfy is given by equation 52. 

Referring to equations 52, 54 or 55, it is seen that the effect of flight to 
flight variations in the PSD of the acoustic pressure jn the predicted com- 
bined load depends (1) on the relative magnitude of P? /v with respect 

p S p (Xj,W) 

to the excitation pressure S p (W) without flight-to-f light variations 
considered and (2) on the magnitude of the transient load considered by 
itself (i.e., if the transient load X(t) is very high compared to the 
acoustically induced load, then consideration of flight-to-flight variations 
in the acoustic load may not oause an appreciable change in the. combined 
load). Similarly, when the aooustically Induced respbnse is based upon the 
M percentile PSD of the acoustic pressure, given by equation 52 
is substituted for o In the equations for calculating First-Passage 
probability. 
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7.0 OOttCLUSIdN 


An analytical technique has been developed for predicting the combined loads 
during STS launoh due to transient meohanioal Inputs and aooustlcally lnduoed 
random vibrations. The prooedure assumes that the individual aooustlo load 
and meohanioal load have been previously obtained by ourrently employed tech- 
niques alluded to herein. The prooedure is based upon three assumptions: 

(1) That the system (structure) is linear thereby permitting superposition of 
the two loads; (2) that the aooustlo load is Gaussian; and (3) that the 
aooustlo load is stationary. During STS launoh the third assumption of the 
acoustic load being stationary is not physically Justifiable, but as the true 
evolutionary nature of the acoustic pressure field may not be known, and in 
fact the acoustic forcing function often is taken as a stationary envelope of 
the evolutionary acoustic SPLs, the assumption of the acoustio load being 
stationary is made as a pragmatic analytical expedient. Should the 
evolutionary nature of the aoouatie pressure be known, this non-stationary 
forcing function could be approximated as being quasi-stationary in temporally 
contiguous intervals. The procedure results in a predicted probability 
distribution function for the combined load; the combined load has been shown 
to be represented as a non-stationary Gaussian random variable with 
time-varying mean, finally, a technique for obtaining combined design loads 
under a first-passage failure criterion has been developed and empirical 
verification of thi proposed technique provided. A summary of the procedural 
steps leading to a First-Passage failure design level are listed in Table VII. 

It is recommended that the following additional work be performed. 

(1) Develop and document computer oOde "run streams" for performing 
first passage failure analyses on Univao, VAX, and CDC machines. 

(2) Empirically verify analytioal techniques for predicting aooustlo 
indueed loads and their resultant strains by employing the data, 
obtained in the experimental investigation, Seotion 4.0. 
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TABLE VII 

A SUMMARY OF PROCEDURAL STEPS LEADING TO A 
FIRST-PASSAGE FAILURE DESIGN LEVEL 


1. Identify acoustic pressure SPL 

a. or SPLs if using quasi-Stationary approximation 

b. if accounting for flight-to-flight variations select desired M 1 
percentile SPL (subjective). 


2. Using the selected SPL (or SPLs) calculate the standard deviations of the 
acoustically induced response, <T y , and the time rate of change of the 
acoustically induced response, <*•, by a random reponse analysis of the 
structure. 

3. Select desired transient and perform responre analysis to obtain 
mechanical transient response, X(t), and time derivative of transient 
response X(t). 

H. Select candidate barrier levels, A. 


5. perform First-Passage failure analysis (equation 37) to obtain the 
probability distribution function of first-passage failure. 

6. Either: 

a. select the desired probability (subjective) to which the structure 
is to be designed and obtain from the probability distribution 
funotion the corresponding design level, or 

b. select the design level and. obtain the corresponding probability. 
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APPENDIX A 


PRELIMINARY CONSIDERATIONS ON 
THE ACOUSTICALLY INDUCED LOAD 


The purpose of this appendix is to develop Some of the salient features of a 
random prooeds (the aooustioaliy induoed load) with which the reader should 
be familiar and that are employed in the text of this document* 



original PAGB 19 

OF POOR QUALITY 


The aooustioally induced load ia taken as a random population, assumed to be 
stationary and Gaussian and of the form: 


V” * SA n SIN< V * V 


or, at a given time, t^ 


W ■ S 4 n StN( Vl * V 


where: 


A n s Amplitude of the n frequency component 

i u 

« n s Circular frequency of the n frequency component 

4 s Phase angle of the n l frequency component for the j member 
of the population 


Each element, Yj(t), of the ensemble has the same power spectral density 
-(PSD) obtained previously from a random response analysis employing the 
acoustio field in the STS cargo bay as a stationary random excitation* The 
elements of the ensemble differ and are deterministically unknown because the 
phase angles differ and are unknown. In this formulation the load due to 
aeoustica is represented as a number of sinusoidal components with random 
phase angles uniformly distributed over the range (-IT , -fir ). As the phase 
angles are assumed to be uniformly .distributed between -tt and +ir , the 
probability density function of is flat at the level, iir, in that 
interval. 


It is known that for some function. of 4 n , say ti(4 n ), the expected value 

of G( 4 ) is defined by: 
n 


lG(0 n )l = f 0(v n’P'V^n 
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where P( 4> n ) ia the probability density function of ^ ri . 

In particular the following will prove to be useful. 

so If 

E|c0S* n | * Jcos* n P( *n )d *n * Jjir* d *n * *# [ slN ^n| n B 0 
• /siN*„ r -in 

* / SIN *n P( *n ,d *n ' j ST - 2 d *n * U l' cos *nJ/ 0 


E SIN*, 


- so 


- It 

v 


I * so 2 

c °S 2 0 n I = jc OS 2 0 n P(0 n )<*0 n a J ^ 0S 0 n S * n |-^L + i SIN 2* n J 

~ * - TT - IT 

= 4^/2 + n l 2 ] • i tA-4) 

E |s>N 2 * n ] = fsm 2 \ «♦„><**„ = f^!^n d*„ '»n[*n -iSlN2* n j 

-* -IT 

= Jirty 2 + */* 1 * * 

so ^ 

= JsiN0 n COS<^ n P(«l> n )d0 n * J|Tt(SlN^ n COS0 n )d^ n 


E SIN0 n COS0 n 


* iir[i SIN 2 ^]" = 0 

• IT 


In a similar fashion, we consider the second order probability density 
funetiort between two phase angles, ^ n » and As before the phase angles 
are assumed to be uniformly distributed between -tt and *v (i.e, -it S 
S it and -tt < < +tt, euch with uniform distribution). Hence the 

second order probability density function between and is a riat 
surface of level ^tt". By definition, two random variables are 
independent if: 


5b 
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(A-5) 


P 2 <*„•*„> * P 1 <*„> P 1 <V 


whioh is the case at hand (i.e. -^"'n 2 * »ir x iir ). Thus, we have 

taken the phase angles to be not only uniformly distributed, but statistically 

Independent* 


It is known that for some funotion of ($ R , 4 m ), say Q($ n , 4> m )t the 


expeoted value of G($ n , $ m ), 1® defined by? 

e ia<v» m )l »V p !'V. )d 'n 4 * 1 


m 


In particular, the following will prove to be useful. 


E 


/ < coso n fCOSO 

~nr **nj-nr* 


d#„ « 0 
m 


E tSIN* n SIN0 m 


/ SIN0 

TT 1 d *n 


r 


SINS 


J 


m 


dO = 0 


m 


E [81N* n COS.J = j d,„ j\ 


COSO 


2 tt. 


— = ® 


(A-6) 


(A-7) 


We had the family of loads due to aooustic excitation to be of the form: 

Y( V 8 sA n SIN ( Vl + *n> (A “ 8) 

where Y(t^) is a random variable because is a random variable. A R 
is obtained from a frequency analysis of response to aooustic exoltation, or, 
by ah identity of trigonometry 



m 
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(A-9) 


(A-10) 


Y(t l ) 8 £A n [ SIN < w n V C0S 4 C0S(W nV S1N *n] 

taking expeoted values 

E « SA n JE jsiNW^COS ♦ E [cOSW^SIN * J} 

and as SIN W t. and COS are uneorrelated 
n 1 *i . . 

and as COS W t, and SIN 6 are uncorrelated 
n X n 

E tY(tj)] - IA„ | E t8IN(W n tj)) E [COS» n l ♦ E tCOSW^)! E lSIN* n l ] «-H> 


but as W t- is a deterministic quantity 
n x 

E [S!N(W n t 1 )J * SIN(W fl t^) 

E ^COS(W n t 1 )j s COS(W n t 1 ) 

and we have already shown that 

E £cOS ^ n J s E £siN 8 

hence 

E jY(t x )j 8 2A n |(SIN W^MO) + (COS W^HO)} a 

M 


E 


8 0 


(A-12) 


(A— 13) 


(A-1M) 


Thus, the expeoted value (mean) of the response. due to aooustio loading is 


aero. 
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Now, squaring the response due to acoustic exoltation gives 
y2(t l ) 8 n A rt ISIN(w n t 1 )COS0 n + COS(W n t i )8IN0 n J 2 

* 2 I I V. ISIN < W n ‘l>COS» n + COS(W n t j)8IN 0 n J 

x[SIN(W m t 1 )COS0 (n + COS(W m t 1 )S!N0 m ] 


or expanding and taking expeoted values 


ECY 2 ^)] s J aJ [(StN^tj) E(C£sV^) + (COS^tj) 
+ 2 SINW^JCO&W^) E(SlK0^OS0 n ) ] 
+2 -n S m A n A ro [siNCw^SINCW^) E(COS^6oS0 m ) 

+ SIN(W n t 1 )COS(W nJ t 1 ) E(COS0^SINp m ) 

+ ElNCW^pcOSCW^p E(SIN0^COS0 m ) 

+ COS(W n t 1 )COS(W m t 1 ) E(SIN 0^SIN 0 m )| 


and E 


|\ 2 (t^)J reduoed to 

E |Y 2 (t 1 ) 1 * £ A n^ SlN^tj +/os 2 W n t J 




S A 


which is Parseval's formula for a Fourier series (i.e. the mean square value 
equals 1/2 the sum of the squares of the component amplitudes)-. 
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In praotice f the meah-aquare aeoustioally induded reaponae, E [l 2 (t)] ia 
the area under the one-alded acoustically induced reaponae PSD, S (W), and 
i© obtained by integrating the aaroe over frequency , i*e# 

E [* 2 (t)J ■ / S y (W)dW (A-18) 

Further, as the expected value (mean) of the acoustically induoed response is 
zero, the mean-square acoustically induoed response equals the varianoe of 
that response.. 


E 




(A-19) 
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EMPIRICAL VERIFICATION OF 
UNDERLYING ASSUMPTIONS 

Having aohieved rather good agreement between the theoretioal First-Passage 
prediction and the empirioal data, it la not aurpriaing that the aaaumptlons 
upon whioh that prediction is baaed were reaaonably valid for the teat con- 
figuration. The prediotion of the oombined response is based upon three 
assumptions: (1) that the structure is linear, thereby permitting superposi- 

tion of the two loads; (2) that the acoustlo load is Gaussian; and (3) that 
the aooustlo load is stationary. Henoe , the oombined response is represented 
as a non-stationary Gaussian random variable with time-varying mean (t>e 
time- varying mean being the response to the transient). 

To oonfim the validity of this representation, empirioal ensemble statistics 
for the oombined loading responses (i.e. statistics across the ensemble of 100 
samples) at various arbitrary time points were compared with the predicted 
ensemble statistics based upon equation 16. 

p*( V = #» t > ♦ •c p * I 

Figures B-l, B-2, and B-3 depict such, comparisons for the aocelerometer #6 
data, speotrum A/transient 1 loading condition of 4-second duration., at the 
times t^ s l second, t^ * 2 seoond, and tj s 3 second, respectively. 

The agreement, if not excellent, is considered adequate. Further, if a 
temporal average of the ensemble statistics over the 4-second duration (i.e., 
over 6000 digitized time points) is made, excellent agreement is found as. 
seen in Figure B-4. Similar results for the other combined loading condi- 
tions for aooelerometer #6 and strain gage #2 are presented in Figures B-5 
through B-20. Taking the oombined response as a non-stationary Gaussian 
random variable with time-varying mean appears to- be a reasonable assumption 
that is essentially oorreot. 
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Oaussl en Distribution Comparison 
Accelerometer 6/Translent 1 ♦ Spectrum 
A/T-2 Second 
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GAUSSIAN DISTRIBUTION COMPARISON 



Gaussian Distribution Comparison 
Accelerometer 6/Transient 1 + Spectrum 
A/Time Averaged 







GAUSSIAN DISTRIBUTION COMPARISON 
ACCELEROMETER 6 / TRANSIENT 2 + SPECTRUM B / T-1 SECOND 



Gaussian Distribution Comparison 
Accelerometer 6/Transient 2 + Spectrum 
B/T-l Second 






GAUSSIAN DISTRIBUTION COMPARISON 
ACCELEROMETER 6 / TRANSIENT 2 + SPECTRUM B / T-2 SECOND 
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Fig. B— 6 Gaussian Distribution Comparison 

Accelerometer 6/Transient 2 + Spectrum 
B/T-2 Second 



GAUSSIAN DISTRIBUTION COMPARISON 
ACCELEROMETER 6 / TRANSIENT 2 + SPECTRUM B / T-3 SECOND 



Gaussian Distribution Comparison 
Accelerometer 6/Transient 2 ♦ Spectrum 
B/T-3 Second 
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GAUSSIAN DISTRIBUTION COMPARISON 
ACCELEROMETER 6 / TRANSIENT 3 + SPECTRUM B / T-.75 SECOND 
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Fig. B-ll Gaussian Distribution Comparison 

Accelerometer 6/Transient 3 + Spectrum 
B/T- .75 Second 
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GAUSSIAN DISTRIBUTION COMPARISON 
STRAIN GAGE 2 / TRANSIENT 1 + SPECTRUM A / T-2 
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Fig« B— 14 Gaussian Distribution Comparison 

Strain Gage 2/Transient 1 ♦ Spectrum 
A/T-2 Second 
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GAUSSIAN DISTRIBUTION COMPARISON 
STRAIN GAGE 2 / TRANSIENT 2 + SPECTRUM B / T-2 SECOND 



Gaussian Distribution Comparison 
Strain Gage 2/Transient 2 + Spectrum 
B/T-2 Second 
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